Cometary Volatiles and the Origin of Comets 


Michael P. A'Hearn 1 , Lori M. Feaga 1 , H, Uwe Keller 2 , Hideyo Kawakita 3 , Donald L. 
Hampton 4 , Jochen Kissel 5 , Kenneth P. Klaasen 6 , Lucy A. McFadden 7 , Karen j, Meech 8 , 
Peter H, Schultz 9 , Jessica M, Sunshine 1 , Peter C. Thomas 10 , Joseph Veverka 10 , Donald 
K. Yeomans 6 , Sebastien Besse 1 , Dennis Bodewits 1 , Tony L. Farnham 1 , Olivier 
Groussin 11 , Michael S, Kelley 1 , Carey M. Lisse 12 , Frederic Merlin 13 , Silvia Protopapa 1 , 

Dennis D. Wellnitz 1 

department of Astronomy, University of Maryland, College Park MD 20742, USA; 
institute for Geophysics and Extraterrestrial Physics, Technische Universitat 
Braunschweig, 38106 Braunschweig, Germany 
department of Physics, Kyoto Sangyo University, Kamigamo JP Kita-ku, Kyoto 603- 

8555 Japan 

4 Geophysical Institute, University of Alaska Fairbanks, 903 Koyukuk Dr, Fairbanks, 

AK, 99775 

5 Max-Planck-Institut for Solar System Research, Max-Planck-Strasse 2, 37191 
Katienburg-Lindau , Germany (retired) 

6 Jet Propulsion Laboratory/Caltech, Pasadena, CA, 91109 
7 NASA Goddard Space Flight Center, Greenbelt MD 20771 
8 lnstitute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822 
department of Geological Sciences, Brown University, Providence, RI 02912 USA 
10 Department of Astronomy, Cornell University, Ithaca, NY, 14853 
11 Laboratoire d'Astrophysique de Marseille, Universite de Provence and CNRS, 

13013 Marseille, France 

12 Space Department, JHU-APL, 11100 Johns Hopkins Road, Laurel, MD 20723 
13 LESIA, Observatoire de Paris, Universite Paris 7, Batiment 17, 5 place Jules Janssen, 

Meudon Principal Cedex 92195, France. 


ABSTRACT 

We describe recent results on the CO/CCh/TUO composition of comets and compare 
these with models of the protoplanetary disk. We argue that the cometary observations 
require reactions on grain surfaces to convert CO to CO 2 and also require formation 
between the CO and CO 2 snow lines. T his then requires very early mixing of 
cometesimals in the protop lanetary disk analogous to the mixing described for the 
asteroid belt by Walsh and Morbidelii (20 1 1 ). We suggest that most comets formed in 
the region of the giant planets, the traditional source of the Oort-cloud comets but not of 
the Jupiter- family comets. 


1. INTRODUCTION 



Our current understanding of cometary dynamics (e.g., Rickman 2010) implies that 
Jupiter-family comets (JFCs) formed primarily beyond Neptune, while long-period 
comets (LPCs), coming to us via the Oort cloud, and likely the source of many Halley- 
type comets (HTCs), formed both in the region of the giant planets followed by direct 
ejection to the Oort cloud or formed anywhere out to beyond Neptune and reached the 
inner Oort cloud via the scattered disk. Recent results on the volatile composition of 
comets, both in situ (A' Hearn et al. 2011) and via remote sensing (Ootsubo ei al. 2011) 
suggest a more complicated scenario in the earliest stages of cometary formation and 
evolution. All comets may have formed in the region of the giant planets. These results 
also call into question the traditional composition deduced for volatiles condensing from 
the gas phase, notably HiO, CO 2 . and CO. in the protoplanetary disk (e.g., Williams and 
Cieza 201 1), although very recent results on the composition of icy grains in quiescent 
clouds suggest a resolution of some of these problems (Noble et al. 201 1). 

Belton and Melosh (2009) have argued that the heterogeneity of the abundances of 
volatiles in JFCs should be explained as a result of evolution in the inner solar system. 
Here we explore the alternative, namely that the heterogeneity is largely a result of the 
original formation process. A similar conclusion has been reached in other recent 
investigations that will be discussed below. 


2. COMETARY VOLATILES 

The relative abundances of HjO, CO 2 . and CO in cometary outgassing have been the 
subject of very limited studies due to the difficulty of the observations, particularly for 
CO 2 . which can only be observed beyond Earth’s atmosphere. The data are summarized 
in Table 1 where we include some dynamical parameters together with the relevant ratios 
of abundances. The first measurements of all three species were for comet 1 P/Halley 
from the in situ mass spectrometers. Those data have had a variety of interpretations but 
the abundances we use here are the ones found by Rubin et al. (201 1 ), who do not require 
an extended source for the CO. Subsequent measurements included an indirect 
measurement via prompt emission in the ultraviolet Cameron bands of CO and direct 
detections with ISO. These measurements have been discussed by Bockelee-Morvan et 
al. (2004), More recent results include measurements from the Deep Impact spacecraft 
of 9P/TempeI 1 (Feaga et al. 2007a, 2007b) and of 103P/Hartley 2 (A’Hearn et al. 2011, 
Feaga et al. 201 1) coupled with measurements of CO in Hartley 2 using HST (Weaver et 
al. 201 la. b). Meanwhile, the Japanese AIC4RI satellite has been used to carry out a 
survey of the three species in about two dozen comets (Ootsubo et al. 2011), although in 
most cases the results for CO are only upper limits, which still provide interesting 
constraints, We do not include here several measurements of only CO w ithout CO 2 . 

Observations of the coma may be misleading and/or confusing due to the presence of 
volatiles, principally 1TO in the present context, in the form of ice rather than gas. as 
seen close to the nucleus in Hartley 2 (A’Hearn et al 2011) and much further from the 
nucleus in the outburst ejecta of l7P/Holmes (Yang et al. 2009). If one observes with a 
sufficiently large aperture that all the ice has evaporated before reaching the limit of the 
field, the measurements of gaseous water should correctly reflect the total rate of release 
of water. A wealth of data suggests that the different volatile ices are not mixed at the 



microscopic level but at some macroscopic level, the scale of which is not known, 
although we note that the volatiles excavated by Deep Impact at the scale of tens of 
meters were in the same relative abundances as those in the ambient outgassing of that 
comet. For this paper, we use the total release of water, both as ice and as gas, as the 
relevant measure for abundance ratios. We consider only observations within 2.5 AU of 
the sun, distances at which ILO-ice is normally expected to entirely sublimate within the 
field of view of typical remote sensing instruments. This is also the distance inside which 
the production rates of OH (from H 2 O) and CO were parallel in comet Haie-Bopp 
(C/1995 01) (Bivere/ al. 1997). The increase of abundances of other volatiles relative to 
water at larger heliocentric distances is easily seen also in the multiple observations of 
individual comets in the complete dataset of Ootsubo et al. (201 1, see also Table 1). 

The abundance ratios for many comets are given in Table I . We have considered 
only comets for which there are measurements or upper limits for CO 2 . This table 
includes measurements at large heliocentric distances. Ootsubo et al. (201 1) have 
explored the variation with distance of observation and shown that the ratio [CO 2 /H 2 O] 
varies over a wide range, w ith a lower limit that is independent of heliocentric distance 
and an upper limit that increases systematically with distance. Here we are interested 
more in long-term evolutionary effects and in Figure 1 we have plotted [CO 2 /H 2 O] versus 
perihelion distance, a parameter correlated with evolution through the maximum heating 
around the orbit. There is no trend among the JFCs and among the HTCs/LPCs, the lone 
object that might suggest a trend is Halley, at 0.58 AU and the only HTC in the sample. 
Note that the range of values for JFCs completely encompasses the range for LPCs. 

Figure 2 shows that ratio for the LPCs + Halley (a much smaller sample), plotted against 
the reciprocal semi-major axis of the orbit on approach, a useful proxy for the number of 
perihelion passages in the inner solar system. Again there is no trend except for Halley 
itself being an outlier with respect to the LPCs. The lack of trends in these diagrams, 
particularly the lack of any pattern for LPCs relative to JFCs, suggests that a primordial 
explanation for the variation should be considered, even though evolutionary effects 
likely play some role. There is no way, with present data, to know whether the low 
values for Halley indicates a trend or is an artifact of smalt-number statistics. 

Meanwhile, the ratio [CO 2 /COJ varies by much more than two orders of magnitude in 
the limited sample for which direct measurements are available. There is not much range 
among the LPCs and HTCs wTiere it is directly measured, but the lower limits for some 
LPCs require that the range be well over one order of magnitude. There is some 
indication from the tower limits that perhaps [CO 2 /CO] might be highest for the 
dynamically new comets. The two JFCs for which both species are detected yield ratios 
more than two orders of magnitude apart and spanning the entire range (including lower 
limits) for LPCs and Halley. This does not seem consistent with depletion of the highest 
volatility species through perihelion passages, nor does it seem consistent with bond- 
breaking by cosmic rays while in the Oort cloud. This also suggests that the ratio may be 
primordial. Mumma and Chamley (201 1) present results from near-infrared 
spectroscopy on [CO/HjOj (their Figure 8) in a dozen comets. Although they do not 
have results for CO 2 , the range of [C'O/TTO] is well over an order of magnitude, with 
JFCs and LPCs overlapping in most of that range. This also suggests similar origins. 



When considering the abundance ratio one must ask whether the relevant ratio is the 
ratio of species leaving the nucleus as gas or the total ratio of species, whether gaseous or 
solid when leaving the nucleus. For some problems, such as studying the non- 
gravitational accelerations or torques on the nucleus, the ratio of gases is relevant, but for 
bulk chemical composition the total amount leaving the nucleus is more relevant. If the 
nucleus is globally heterogeneous and the different volatiles are coming from 
fundamentally different areas within the nucleus, the measured abundances may average 
over different cometesimals. The fact that the excavation by Deep Impact at Tempel 1 
led to volatile abundances in essentially the same proportion as in the ambient outgassing 
(A' Hearn ei al. 2008) argues that at least to depths of tens of meters the ices are mixed 
(macroscopically, not microscopically). These results are very hard to explain in an 
evolutionary scenario and we consider the primordial scenario more likely and therefore 
worth studying in more detail. 

3. PRIMORDIAL. ABUNDANCES AND THE PROTOPLANETARY DISK 

We first note that most traditional models of the protoplanetary disk predict 
abundance ratios CO 2 /CO « 1 (typically as low as 10' 4 ) reflecting gaseous relative 
abundances and direct condensation from the gas. However, recognition of much higher 
ratios in the ices of interstellar dark clouds has led Garrod and Pauly (20 1 1 ) to 
theoretically consider reactions of OH with CO on the surface of icy grains that lead to 
comparable abundances of CO and CO 2 in interstellar dark c louds. Noble el al. (20 1 1 ) 
have obtained similar results in the laboratory. We conclude that these surface reactions 
are critical to the format ion of comparable abundances of CO and CO 2 in cometary 
nuclei. 

We further suggest that the wide range of abundance ratios in comets implies that the 
comets formed generally between the CO 2 sand CO snow lines in the disk. The location 
of these snow lines differs among models and even varies with the evolution of the disk 
for any single model. For some recent models of the protoplanetary disk, e.g . . that by 
Dodson-Robertson el al. (2009), these tw'o snow lines are in the region of the giant 
planets. It seems unlikely that inclusion of grain-surface reactions would have major 
effects on the thermal profile of the disk. Thus we suggest that most comets formed in 
the region of the giant planets, the traditional origin for Oort-cloud comets but not for 
JFCs. 

It is widely understood that the giant planets migrated during the early stages of the 
solar system although ideas on the details of that migration are changing rapidly. 

Recently Walsh and Morbidelli (201 1), in order to explain the distribution of asteroidal 
orbital properties, have invoked migration of Jupiter and Saturn beginning much earlier 
and much closer to their present-day orbital positions than in recently preceding 
dynamical models. This migration mixed asteroids both outward and inward as the 
direction of migration changed and ultimately mixed up most of the asteroids in the outer 
portions of the asteroid belt. We suggest that a similar process must have occurred for 
the comets, although the role of specific planets is not clear without more work. Some 
comets were ejected directly to the Oort cloud as classically understood, but many others 
must have been ejected to the scattered disk and thus provided a source for the JFCs. The 
details of this process require a detailed, quantitative simulation that is beyond the scope 



of this paper, one involving grain reaction rates as well as specific planetary migration 
rates since the time scale is short. 

We note that such mixing could also explain heterogeneity within individual 
cometary nuclei since this very early planetary migration could have occurred while the 
eometesimais were still accreting into nuclei. A’Hearn el al. (201 1 ) and Feaga el al. 

(201 1) have discussed the variations in relative abundances in Hartley 2 as it rotates and 
this strongly suggests that the larger lobe of the nucleus has much less volatile content 
than the small lobe. It is not yet clear how the relative abundances of the three ices varies 
because a detailed, quantitative model of the sublimation of grains in the coma is 
required. Because this comet is in excited state rotation now and the observed changes in 
rotational state during the last apparition imply different rotational states earlier in its 
history, seasonal effects should be minimal, suggesting that the two lobes of the nucleus 
likely formed at different heliocentric distances. The obvious spatial heterogeneity in the 
outgassing of Tempel 1 (Feaga et al. 2007a), on the other hand, could be either a seasonal 
evolutionary' effect or a primordial effect. A variety of other earlier measurements by 
remote sensing have reported chemical changes with rotation suggesting that other 
comets also have nuclear heterogeneity, although in those cases it is not possible to assess 
the relative importance of primordial abundances, seasonal effects, and long-term 
evolutionary effects. 

Mixing of this sort also appears to be needed to explain the D/H ratio of comet 
Hartley 2 as pointed out by Hartogh el al. (201 1 ). They measured the ratio to be equal to 
SMOW (Standard Mean Ocean Water), i.e.„ a factor two less than in the few 
LPCs+Halley for which it is measured. This suggests that at least Fiartley 2 was formed 
in a region with D/H similar to that on Earth. If JFCs formed in the Kuiper belt and LPCs 
formed in the region of the giant planets, the differences should go in the opposite 
direction according to Hartogh et al. Furthermore, Hogerheijde et al. (20 1 1 ) have 
argued that the ortho-para ratio of water in a protoplanelary disk, which is much lower 
than that in comets in our own solar system, also requires comets to contain material that 
condensed in a wide variety of locations in the protoplanetary disk. 

4. SUMMARY 

Observations of the abundances of the volatile species H 2 O. CO 2 . and CO in comets 
show a wide range of abundance ratios in both JFCs and LPCs/HTCs. This seems 
unlikely to be the result of evolution through successive perihelion passages. We suggest 
that both classes of comets, or at least the eometesimais that might later have assembled, 
formed between the COi and CO snow lines, thus allowing the wide range of abundance 
ratios. They were subsequently scattered by planetary migration to lead to the present 
dynamical types of comets, much as the asteroids were mixed up in the dynamical model 
of Walsh and Morbidelli (201 1 ). This could also explain variations within individual 
nuclei since this migration takes place earlier than previously modeled and is relatively 
fast. 
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FIGURE CAPTIONS 


Figure 1. Relative Production of CO 2 as a function of perihelion distance - 
logarithmic scales. jFCs show no trend with perihelion distance and span more than 
the entire range of relative production rates by LPCs. The point in the lower left of 
the diagram is lP/Haliey, the only HTC in the sample. 


Figure 2. Relative production of CCF? as a function of original reciprocal semi-major 
axis for LPCs and HTCs - logarithmic scales. 1 P/Halley is the unusually low point 
and is the only HTC in the sample. 



Table 1. Abundance Ratios of Volatiles 
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Note: Columns: 1. Object name and ID; 2. Heliocentric distance at time of observation; 3. Perihelion distance at observed 
apparition; 4, Reciprocal semi-major axis on approach to planetary region; 5-10, Ratios of production rates with errors or 
limits as appropriate; 11. Reference for data. 

Abbreviated Names: 1. S-W 1 = Schwassmann-Wachmann 1; 2. Chur-Ger = Churyumov-Gerasimenko; 3. S-L 4 = Shoemaker- 


Levy 4 
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